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PREFACE 

The progress of resesrch on "Three Dimensional Flow Field Inside a 
Compressor Rotor Blade Pessage, Including Blade Boundary Layers" (NASA 
Grant NSG 3266) for the six-month period ending December 31, 1981, Is 
briefly reported here. Dr. J. M. Galmes, Research Associate in Aerospace 
Engineering, has the responsibility for turbulence modelling. Mr, Pouagare, 
a doctoral candidate, hafe assumed responsibility for the development of the 
coiqiuter program to predliit the blade boundary layer. The measurements 
reported were carried out by M. Pouagare and K. N. Sachldananda Murthy. 

B. Lakshmlnarayana 
Principal Investigator 


TABLE or CONTEMTS 

P«f 

NCMENCLATURE Iv 

1. NUMERICAL ANALYSIS OF BUDE AND HUB WALL BOUNDARY UYERS . . . 1 

Introduction . . 1 

Momontun Equation* ... 2 

Nuaarlcal Solution of th« Three Monentua Equation* ...... 4 

Pr***ure Correction Equations ................ 6 

Statu* of the Goaputer PrOgraa ................ 9 

2. TURBULENCE MODELLING FOR BOTH HIGH AMD LOW REYNOLDS NUMBER 

FLOW SUBJECTED TO ROTATION .................. 10 

Introduction . . . . . . . . . ...... 10 

Turbulence Equation* . ............ 13 

Qualitative Analyals of k and e Equation* .......... 16 

k>e Model for High Reynold* Number Flow* ........... 19 

k-e Model for Low Reynold* Nuaber Flow* ........... 22 

Eddy Vl*co*lty Law .... ............. 27 

Conclualons ............... .... 28 

3. MEASUREMENTS OF THREE DIMENSIONAL FLOW FIELD INSIDE A 

COMPRESSOR ROTOR PASSAGE . 30 

Experimental Program 30 

Typical Reault* 30 

4. PAPERS, THESES, AND REPORTS PUBLISHED DURING JULY 1981 - 

DECQIBER 1981 33 

REFERENCES 34 

APPENDIX 1 Transformed Momentum Equations and Jacobian Matrices . 39 

APPENDIX 2 Literature Survey on the Effect* of Curvature and 

Rotation on Turbulence Structure ............... 48 


FIGURES 


61 


Iv 


WOMEKCLATURE 


♦ 


constants 


^63 

Modeling constants 

fc<*r *i«> 

dunping function for low Reynolds number flows 

«lk 


awtrlc tensor 

1» J 


Indices In the C»n»R coordinates* respectively 

J 


Jacobian of the tranafonutlon 

^ - i 8^“ ^ 

turbulent kinetic energy 

1 


length scale 

p 


pressure 

p 


mean pressure 

p' 


fluctuating pressure 

PS 


pressure surface 

Ps 


static pressure 



total pressure of the relative flow 

**TOTAL 


Vl < 

**STAT1C 


< 

r,0,z 


cylindrical coordinate system 

R 


radial distance normalized by the tip radius 

Re 

u'^ 


Reynolds number 

Re ■ “v" 

D 

*^T " ve 

» 

local Reynolds numbers 

» 

e._ 




generalized gradient Richardson number 
3R/3r 

"t.j 

R 

r 


s 

ss 


s 


IJ 


u 


I 

> 


1 


t 




“i“j 
u,v,w 


V 

w 

WR 

WT 

WZ 

Y 




e - 2V 


ipq 


n 

r\c 


blada ipaoing 
sue Cion surf act 
sCraln Cansor 
fluctuating valoclty 
time 

Man contravarlant valoclty 
Reynolds tensor 
blade tip speed 

velocities In radial > tangential, and axial direction, 
respectively 

Kolmogorov velocity scale 
relative velocity 
radial velocity, Vf|./U^ 
relative tangential velocity, Wq/U^ 
axial velocity, W /U^ 

contravarlant and covariant coordinates variables 

tangential distance measured from camber line (y - 0 on 
suction side, Y - S on pressure side for the data In the 
passage for Che wake, Y ■ 0 Is the center of the wake) 

axial distance normalized by the blade chord (Z ■ 0 at 
the leading edge) 

constants 

Kronecker tensor 

turbulent dissipation rate 

permutation tensor 

Kolmogorov length scale 

an/rse 


vl 


U 

% 

5 f n » R 


h 

p 

n 


■olccular vlicos.Uy 
turbulent viscosity 

trsnsfornsd coordlnstss In the strssiMlss, nornsl, and 
radial directions, respectively 

35 / 3 * 

35/r36 

density 

Modeling constants 
angular velocity 

contraverlsnt component of angular velocity 


Subscripts 

l,j,k,l,Bi,n 

(p) 

r 

s 

T 

z 


Indices 

value at the wall 
radial component 
static 
total 


e 


axial component 
tangential component 


1. NUMERICAL ANALYSIS OF BLADE AND HUB WALL BOUNDARY UYERS 


Introduction 

The flow In • turboMchlncry blade paaaa^e hae a predominant flow 
direction. The vlecous dlffuelon In the atreanwlae direction le uaually 
email and the elliptic Influence la tranamltted upetream through the 
preaeure field. So the atreanwlae dlffuelon teima In the Navler«>Stokee equa» 
tlone C|tn be neglected In comparlaon to the diffusion In the radial end 
tangential directions. Starting with a guessed pressure field, It Is 
possible to converge on the full elliptic solution by Iterating between 
a parabolic solution and an Iteration of the pressure field, ^ The main steps 
of the calculation are given below. It Is assumed that the properties at 
the streamwlse station 1 are known and solution at the streamwlse station 
1 <f 1 Is sought. 

1. The three momentum equations are solved simultaneously to get 


U, V, HujU J 




2. The energy equation and the perfect gas relation are solved to get 
the temperature T and the density p . 

3. The streamwlse pressure gradient Is corrected by the use of mass 
conservation across the cross-section. 

4. The cross-flow velocity components, determined from the momentum 
equations, are corrected by Including In them the irrotational velocity 
corrections that satisfy the continuity equation. This step requires the 
solution of a polsson equation for 4 defined later. 

5. A polsson equation Is solved to determine the transverse variation 
of pressure . The equation for the pressure is derived by appropriately 
differentiating and sunnnlng the two transverse momentum equations. 
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6. Staps 1 to 5 art rapaatad until convarganca la achlavad. 

7. TUrbulanca cloaura aquatlona ara aolvad to gat a naw valua for 
tha affactlva vlacoalty. 

8. Solution la Barchad to tha naxt atraamwlsa atatlon* Tha flow chart 
of tha conputar program la ahown on tha naxt paga. 


Tha Momantum Equations 

1 

The Navlar-Stokaa aquatlona In cylindrical coordlnataa arc tranaformcd 
Into the computational domain neglecting vlacoua diffusion tezma In the 
direction 

The transformed equations can be written In the form, 


JE . in 4 . 3G _l[> . ia X el 

35 ^ 3rt ^ 3 r ^ ^ Re|_3R ^ 3n * J 
The derivation of this equation Is given In Appendix 1. 

nn 

The column vector of unknowns Is g |V 

W 


( 1 ) 


The column vectors E, F, G, C, P, Q, S are given by 
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Flow Chart for the Computer Program 




f; 


4 


1 


i: 


PUW/R 

2PVW/R + 2pVM 

L »■ - PJ^^R - 2pm 


J 




% 

i n » is lE fn ^ X « ^ W> 

3 U Rr 3^ - -f [r\Q -gjj' n, 



Pe 



n 

^2 

: J 


f'^e 

"•1 


r 

^2 

j J 
2 

1 1 I 

Pe. 


H 

i 

J J 


S - 


E fo M X n 3Wl 

R L^r 3 r ■•' 3nJ 

iEfn is . « is.v] 

r L^e 3n + 3 r " rJ 

^Fr iS- n i3^ Wl 
R L r 3R 3i rJ 


+ 2p^ii 

^ J R 


Numerical Solu tion of the Three Momtntum Eauatlona 

The numerical scheme is based on the Linearized Block Implicit Method 
of Briley and McDonald (1979). Equation 1 Is written In the form 


- D(g) + S(g) 


( 2 ) 


s 


wh«r« H(g) - E(|) 

D(.) - J.rR iS + ial .iE.ic 
R«lff ary an aR 

s(f) - S(g) - C(g) 

Thft flnlt* difftrtnca fom of oqiuitlon (2) it 


(H 


1+1 


- nh/U - e[i(g) + s(g)J^‘*‘^ + (I - ©|i)(g) + t(g)]^ 


vhoro 0 Is ths wslghting factor. 

Tht llnssrisstlon is dons at follows: 


- H* + (f)‘ A5 . H* . ^)‘ If A5 - h‘ . - ,‘) 

s‘« - . (I)" A5 . S* + (!)*(,*« - ,‘) 

. + (||j‘ A5 . D* 


* f|2)'(.‘« - .‘) 


- D' + 1-^1 - D- + l—l [g 

The llnearizad achaaa is than given by. 


Defining A 


[h] 


- m 


L - -p 


3D 

dg 


and 


. 1+1 Ifl 1 
Ag -g - g 


Equation 3 can be written as 


or 


(A + ASDAg^"^^ i. A5tD + S]^ 


(A + AC(L^ + Lj^))Ag^“*'^ - ACtD + S]^ 


( 3 ) 
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Faetorlclnt th« tchcM band on tha Douslaa-Cunn apllt ona gata 

[A + A5L^)iig* - A?tD + S]^ 
lA Aa JAg^'*’^ - AAg* 

IT .« • 

Tha Jacobian oatrlcaa appearing In aquation 3 ara given In Appendix !• 

Praaaura Correction Eouatlona 

The preaaure la conaldered to be coapoaed of three parte. 

p(5.n.R) - Pi(C,n,iO + Py CO + P^(n,R) (3) 

P^ la the Initial gueaaed preeeure field. Py(^) and P^(n»R) are one- 
and twoodlaenelonal preaaure correctional reapectlvely . 

The velocity coaponenta are decoapoaed aa follova 

u ■ u* + u‘ 

V » V* + V’ (A) 

w - w* + W’ 

U*| V*| U* are the velocity coaponenta which are calculated from the 
moaentum equations. U*. V% W* are the corrections needed for the continuity 
equation to be satisfied. 

The U-moaentuB equation la written as follows. 
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pu*(0* + U*) + -y p(U* + U’)V*| + pu*(u* + U’) + •— p(U* + U')V*j 

, OMHI . ,A[k 

- ^(t »1 + v] +-c{4[7["K + «■> + ^]] 

+ + »■> + ^ t("6 ^«* + “’>) 


The above aquation can bo rearranged to glve» 

5, 


^fr ^ ^ o'""’) + »(^ o«*o’) + ^ 

f 

- - - iFf %) * i :\0 " i I R * 

m 


R r 3R 


Writing the above equation in finite difference form and dropping the 
mixed derivatives and the off diagonal terms we get, 

5. n ^ fK. 


~r ~ pu* + — ^ pv*l u* + u* p 1 

/^5L J ^ ^ j P'' Ji+ij,fc i+lj.k I R Ji+ij,ic m j cj 

»'_l ■' Jl+l,j,k+l ' J 'l+l,j,k-l I " 'l+l.J,k 

+fefng + 1.333 nj]l + ft(n? + 1.333 nfll 

yi 0 *wi+i,j+i.k yi 0 ‘iJi+i,3-i,k 

‘Dl+l,3,k 


- ^lif'e + 1-000 "f 


or 


3Ri 


U' “ A -g^ where A is a known quantity 


(5) 


a 


Using V ■oBsntua squstlon ons gsts 

3p 

V* B whern B is s knoyn quantity 


From continuity aquation, wa gat 


or 


// (u + n v) dR dn ■ » ■ constant 

nR * 


// -^:;:5-(u* + A B -^) dRdn - 6 


nR^, 


X " 35- 


m - // - ^g = =- (U* + n,V*)dRdri 
3P’ ^ nR /l+n6 

// (A + B)dRdn 

nR 


( 6 ) 


( 7 ) 


To correct the cross-flow velocity component the following velocity 
potential is Introduced 


^ ^ p W 

3R J ^ 


( 8 ) 




The continuity equation in the transformed plane becomes 


( 9 ) 


PW ^ 
R 3R 


(tp") + p’ * t p“] - ^fr p’ + T p”) ■ 


( 10 ) 
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fi(W* + W) + p (W. + w)| + 0 (v* + V) + -^ 

+ -^^p(v* + V) ♦•y p(u* + U')] - 0 


p(u* + 0')] 


« «’ ♦ w(t '»'•] + + T o"’) * wfr »''' + ^ PU-) - -i* (12) 

vhcrc A* can bt avalu^tad from V*» W*, p. A* abould ba saro. If continuity 


la aatlaflad. 


Ualng aquatlona 5, 6, 8, and 9 In 12 ona gata 


R 3R ^ 3n ■ ■" ‘ 3cl d " 95 ■*' J ^ ^ 9€ 


In order to darlva the Polaaon aquation for the praaaurai the n-nonantum 
aquation la dlferautlatad with raapact to r\ and the R-nomentum aquation with 
raapect to R> Than they are added together. The reault la 


3 p 




where C contains known quantities. 


Status of the Conputer Program 


The part of the computer program that aolvea the three momentum equation 
has been completed. Coding of the polaaon equations la under development. 
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2. TUltBULENCE MODELING OF BOTH HIGH KEYHOLDS HUMBER 
AND LOW EEtWOLDS HIIHBER PLOWS SUBJECTED TO EOTATtOW 

Introduction 

i 

Th« throu dlnonalonal viscid and turbulant affacts In turbonachlnary 
arc nalnly coucad by the threa dlnenslonal boundary layer on blades and wakes, 
annulus well and hub wall boundary layarc» shock-boundary layer Interaction » 
and secondary flow. Theaa viscous and turbulenca effects play e dosilnant role 
In the study of Improved design, better efficlencyi off design perfomance, 
etc. 

Vfhlla appreciable anount of wotk has been done In understanding the 
Invlacld effects In turboouichlnecy, there has been no Investigation related 
to the blade boundary layers In coepressolr rotora, Because of the coapllcated 
nature of the problen, conpreasor, turbine and rotor fan blade boundary 
layers still continue to be one of the least understood phenomena In turbo- 
machinery. The boundary layers are three dimensional with laminar • transi- 
tional, turbulent and separation zones. The flow field is a function of 
several parameters such as Incidence, solidity, blade geometry, hub/tlp ratio, 
camber, radial and chordwtse pressure gradients, inlet turbulence, Reynolds 
number, Rossby number, Mach number, etc. 

A knowledge of the boundary layer characteristics, both mean and 
turbulence properties, is essential in the prediction of flow behavior in 
these blade passages. The Penn State group, presently involved in the study 
of the endwall flows, has initiated recently a detailed study, on both the 
experimental and analytical point of view, of the blade boundary layers • The 
experimental survey of the flow is in preparation and will be atarted very 
soon. The existing and well proven techniques of single hot wire and x-wlre 
are to be used in measuring the development of boundary layers pn turbomachinery 
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rotor blades. Both MsaurMMinta of Man and turbulant flov flald nro to bo 
parforMd. 

Analysis of throe dlmenslonil boundary layer la equally complex due to 
the fact that additional effects such as Coriolis and Centrifugal forces 
change the structure of turbulence, thus invalidating most of the turbulence 
models that are presently used In computing turbulent boundary layers. In 
the case of turbomachlnery rotor flows, the turbulence Is affected mainly by 
the curvatures of both the blade surface and the streamlines, by the body 
rotation, end In the case of the boundary layer, by the low Reynolds number 
effects. All these effects make the flow to be highly non-lsotroplc . Nore-* 
over, as it can be seen In the experimental results of Castro and Bredshaw 
(1976) and, Johnston, Halleen and lezlus (1972), the curvature and/or the 
rotation may affect the stability of the boundary layer and an augmentation 
or suppression of the turbulence may result. With such phenomena occurring, 

It appears that assumptions based on the well known Isotropic eddy viscosity 
concept should fall badly, as the Reynolds stress tensor is not aligned with 
the mean strain tensor when additional production or destruction of turbulence 
are coupled with the production due to the shear. 

Most of the present models are valid for non- rotating systems. These 
include Jones and Launder (1972) for k-e, Launder, Reece, and Rodl (1975) 
and Lumley and Khajeh-Nourl (1974) for the full Reynolds stress model. The 
effects of curvature in high Reynolds number flows have not yet been accounted 
for properly In the turbulence modeling. We may mention the attempt by 
Gibson and Rodl (1981). As far as the rotation effects are concerned, very 
few attempts have been made to Introduce such effects In the calculation 
schemes. However, In the available turbulence models, the effect of rotation 
is not properly modeled In the transport equations of Reynolds stress or the 
dissipation rate. Moreover, the boundary layers flows are of the low 
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Rtynolds' nunbcr type. If we are intereeted in Modeling the turbulence 
behavior very i>*ar the wall» then it ie neeeaaary to introduce thia effect 
in the Modeling aaeunptione. 

The formulation of the cloture problem will vary depending upon the 
Infomatlon and accuracy deaired. In fact, the calculation proceduree for 
three-dimenelonel vlecoue flows require large Memory storage and large 
computer time to solve the three monentum and continuity equations. 

Moreover, such codes are still in a phase of development and are generally 
tested only for simple cases (e.g., laminar flows). Hence, the introduction 
of turbulence models in these codes necessitate careful attention and simple 
models may permit control of the stability of the numerical scheme. However, 
in rotating turbulent flows, zero equation models fall to represent even the 
gross properties, due to the fact that the length and tine scales are assumed 
empirically. The two-equation model is a compromise between a full Heynolds 
stress model which needs the resolution of seven more transport equations 
as complicated as the Navler-Stokes equations themselves, and the empirical 
models. Therefore, as a first step, attempts are being made to Include the 
effect of rotation in the k-e model. Modeling of the rotation effect and 
the low Reynolds number effect in the turbulence closure equations are 
described in this report. 

A literature survey on both the analytical and the experimental work 
is given in Appendix 2. Some of the important conclusions of this survey 
are: 

1. Only few calculations are available for the calculation of the 
three-dimensional boundary' layer in rotating frames. 

2. No complete Reynolds stress model is available for rotating turbulent 
flows. Very few attempts have been made to account for the rotation effects 
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In th« k-e but none nr« really based on a logical analysis. In fact, 

the Boat up to data results In Bodsllng are those of Raj (1975) and Hah and 
Lakshalnarayana (1980) . The aajor effort, then, should be given to the 
analysis of the dissipation rats equation and to the Reynolds stress equations. 

3. Detailed ■sasursnsnts providing Inforaatlons on the effects of 
both the Rossby and the Richardson nundter on turbulence are not very nuaeroue. 
So, every new result would be of great Interest, particularly If the rotation 
effect can be Isolated froa the other effects. 

Turbulence Equations 

The turboaachinery boundary layers and wake aust be represented In a 
relative rotating fraoie of reference, which Includes both the curvature and 
rotation teras, to ellalnate the effect of periodic unsteadiness. The trans- 
fomatlon of the turbulence equations as well as the aoaentua equations from 
a stationary coordinate system to a rotating coordinate system Is quite 
complicated. Moreover, the necessity for the turbulent model to be as 
general as possible. Indicates that the generalized tensor formulation Is 
the most adequate representation of the set of equations. The equations of 
the mean and turbulent quantities representing an Incompressible flow are 
derived in a rotating frame and are presented In conservative form. 

Continuity equation 

U^i- O u'J - 0 (15) 

Momentum equations 

(p0^)+(pU^uJ),j +2eipqPnPu‘> +p[(fi^ - -(p^+puju'^ - pJ),j 

(16) 


with • 2 ms| for Hcwtoiilan fluid and 


Raynplds atraaa aquatlona 

(puJ3J)+(pl^ bJ),j - -(r3J«J + + puju;;u'i - IljpJ 

♦ «i - \i 

- “k"’’ ♦ ^pq 


(17) 


Klnatlc anergy aquation 


1 »“ “1=1 ■ = 


(pE)-f(pi( U^),j • “(p*u^ 6^^ + pku’^ - u|P^^),J - pu’^u'^ - P^^u|^j 

whara -pu’^u'^ ^1 J " P^^u^ j ■ pe ■ 


(18) 


Dissipation aqtiatlon 

(pe)+(peu^),j • -4y(Sjj^u'-* S^J)-(pcu^) + g"^p e,^j - 4p j 

-4m - 4p - 4pv,"J Si}*' SJ^ „ 

-4V .•^"p.’j^ - 8p Cjp/ S •*''»!’ (19) 

♦ 

The curvature terms are Impllcltely included in aquations IS through 19 . The 
rotation affect appears both explicltely and impllcitely. We are now focusing 
our analysis on the rotation effect principally. 

In the momentum equatloup the rotation appears through the Coriolis and 
Centrifugal forces, but it also affects Impllcitely the results through the 
Reynolds stresses. In fact, looking at equation 17, these Reynolds stresses 


IS 


•r* ixplleltljr affacttd by th« Corlolia forcati, but arc alao lopllcitly 
affaetad by tha rotation through tha triple velocity correlatlona, the 
preaaure velocity correlatlona, the preaaure atraln correlation, the produc- 
tion by tha atreaiea themaelvea, and the dlaalpatlon. 

Following theae renarka It la evident that a full Reynolda atreaa a»del 
ahould give better reaulta than the andels baaed on an laotroplc eddy vlacoalty 
hypotheala> particularly In three-dlnenalonal flowa. 

However, the difficulty to handle aolutlon procedurea for three dlmen- 
alonal turbulent flowa and the particular problea of modeling the dlaalpatlon 
equation lead ua. In a flrat atep, to direct our efforta In developing a k-e 
mdel for both high and low Raynolda nuaAer flowa. Thla andel neceaaltatea 
the reaolutlon of two tranaport equatlona for k and e, and the Reynolda 
atreaaea are related to the taean atraln through an eddy vlacoalty. The 
preaent k-e model cannot account for the anlaotropy of the turbulence which 
exlata in the boundary layer around a blade. A modification to the relation 
giving the eddy vlacoalty la preeently under development. In order to Include 
the effecta of the rotation. 

The rotation alao appeara explicitly In the equation for the dlaalpatlon 
rate e. It aeems [Hanjallc and Launder (1978)] that aome of the major problema 
in calculatlona come from the modeled equation for the dlaalpatlon rate e . 
Hence, there la a need for better analysis of this equation. Following 
Tennekes and Lumley (1972), we may develop a qualitative analysis of the 
equations k and e in order to provide the magnitude of each term and then 
to derive a model for high Reynolds number and low Reynolda number flows. 
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Qualltatlv An>lyl« of k and e Equatioiw 

Equation 18 npraaants tha avolutlon of tha klnatlc anargy k. It la 
Intaraatlng to nota that tha rotation tara vanlahaa Idantlcally, and thla 
aquation la Indepandant of tha coordlnata ayataB. It la claar» howavar, that 
avan though tha Corlolla forcaa hava no dlract affacta In thla aquation, tha 
rotation affacta tha klnatlc anargy principally through tha production by the 
naan atralna and alao through tha dlaalpatlon tara. Tha production tarn 
dapanda on how tha Raynolda atraaaaa are rapraaantad and tha dlaalpatlon tara 
la daacrlbed by aquation 19, which la the exact fora of tha aquation of evolu- 
tion for tha dlaalpatlon rata. Dlffarant authora hava pointed out tha 
difficulty In aolvlng thla equation, and one of the aoat aultabla aathods 
to almpllfy tha aquation la to uae tha high Raynolda nuaber approxlaatlon 
With this asBuaptlon, teraa which are doainant at low Reynolds nuaber are 
ellalnated. At this point, we aay note that equation 19 la not Independent 
of the transforsatlon froa a non-lnertlal fraise of reference to a rotating 
fraae. 

In Tennekes and Lualay (1972), It is shown how ordars of magnitude aay 

be assigned to various correlations appearing In these equations. Instantaneous 

quantities appearing in the correlations are of two types, belonging either 

to the energy containing range of eddies, or to the dissipation range. The 

former has characteristic frequency u'/& (where u' is a fluctuating velocity 

scale, while i Is a length scale of the gross structures) . The dissipation 

rate is of order u* /i. The latter has characteristic frequency v/q (where 

V and q stand for the Kolmogorov velocity and length scales) , with vAl 
1/2 

Re u'/il (with Re u^l/v). The correlation coefficient between two quantities 
from the same range may usually be taken as unity, but the coefficient 
between two quantities from different range is of order of the time scale 


*1/2 

ratio n/v/A/u* ita . In addition* wa nay aaka uaa of tha «ora faalllar 
fact that darlvatlvaa which ara axtamal to corrajLatlona corraapond to acalaa 
In tha anargy containing ranga* whlla darlvatlvaa within tha corralatlon 
corraapond to dlaalpatlon acalaa. Wo wlah to apply thla aort of raaaonlng 
to tha aquatlona for tha klnatlc anargy k and Ita dlaalpatlon rata c. Appllad 
to the dlaalpatlon aquation It la particularly uaaful* bacauaa tha dynanlca 
of thaaa quantltlaa la doalnatad by tha aaall acalaa* and Intaracta only vaakly 
with tha anargy containing addlaa. Procaadlng In thla way* wa may provlda an 
ordar of magnitude for each term of aquatlona 18 and 19. 


Klnatlc Enarav Equation (aquation 18) 

*- - -i 

Tha convective tacma (pE)4’(pk U^)** are of tha aama order aa p 

J <i 


SsdI 


(u*)^ 


~ ' i " j _ 

The production term (^ti* u'*' U. .) la of order p 

i*j * 

(u O ^ 

The dlaalpatlon term (-pe) la of order p — 

In the "dlffualon" term -(p'uj^ 6^^ + pku'^ - uJP^) , the two flrat 


terma are of order p 


(u’)' 


, while the laat one* repreaentlng the dlffualon 


by the molecular proceaaea* la of order p 
high Reynolds number flows. 


(u*)^ ,rl 


Re and la negligible for 


Dlaalpatlon Equation (equation 19) 


The 


(ul)l 


*1 (u* 

convective terms (pe)+(pe U'*) , . are of the same order -p 

J 9.^ 


,J slk 

i ’ II c 


The terms Involving derivatives of mean strain ®ik“’ Is of 

the order p ^ ' ■ Re^ . 

r 

The terms Involving mean derivatives and rotation 




IS 


'■ IV ' JL J 

FoUowlnt Lualty (1970), w« Miy aodcl Ch« quantltl** S'* ul 4 and S' u'r 

* *tj »•* 

and than wa writia 


s'**SijJ s + S/f s^^ j 

and alnca eha naan flow la conaldarad IncoBpraaalbla tha only tarns which 

contrlhuta art tha cacond tans In tha last two ralatlona. Tharafora, tha 

(u')^ -1/2 

first tarn Involvlnt naan darivativas is of tha ordar p Ra and tha 

sacond tarn is of ordar p ^^ 2 ' ^ Ra^^^ (1 + 2 . 

Tha tarns involvlnR tha ttlpla corralation batwaan tha fluctuating 

iL i 

darivativas >4u S' u! .u!' . Following ones again tunlay (1970) it is of 

*sJ 

tha ordar 

(S + 0 ( 81 ^ 8 ^*^ |)) 

1/2 -1 

Tharaf ora It is of ordar p — Re' (B + Ra ) . 

Tha tarns Involving tha corralation with prassura fluctuations 


Wa nay write a poisson equation for this tenn and then the solution of the 
aquation is as follows: [Chou (19A5)] 

dV, 


V s 
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JL fff ri . t ® 

2lt >>> L(l) <1> J 
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^ 4ir |x 3nQ) 
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- PJ 


(D.k" 


mu 31/x 
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dS, 


3n(l)r“(l) 
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This solution Introducoo hlthtr lovol of corrolstion and than, at thla 
laval of closura asauaption va nay naglact thia tarn. 

Tha tama involving tha corralatlon batwaan tha darlvativas of tha 
fluctuating strain 

-4MVg"^ (a +0(8^118^*^ |>) 


Tharafora it is of ordar 


Ra^^ (A + Ra^) 




Th« •’diffusion” tons -(peU’"’)| . *f g Me 9^4* Th« first tsm is of 

A j “j . A 


ordsr P 


(u») 




whlla tha last ona is of ordar p 


/ 


This ordar of aagnltuda analysis is suaMriaad in Tabla 1. 


k-e Modal for Hiah Raynolda Nunbar Flows 

High Raynolda nunbar flows occur ganarally far fron walla, in that 
Casa tha viscous diffusion nay ba naglactad in aquation 18, than only the 
diffusion tarn (pku*^ -f p'u*^),j naads to ba nodalad. For this tarn wa follow 
the proposition of Jonas-Launder and we write (in generalized tensor notation) 

The production tern is defined through the law of eddy viscosity. 

In regard to the dlasipatlon equation, we nay neglect all the terns 

(u’)A 

of order of nagnltude less than p — 5^ ; therefore, the only terns left are 

r 

the convective terns, the "turbulent diffusion" and the production and 
dissipation by the correlations of fluctuating derivatives (Table 1) . These 
two terns have been derived by Lunley *and Khajeh-Nourl (1974) . At this point 



of nodollng* assuaptlons to got thooe two tcxao oro widely based on physicel 
end dlaenslonel anelyeis. In feet, for high Reynolds number flows, the 
dlssipetlon proceiees occur nlmost only et high wave numbers (e.g., the 
smallest structures); therefore, an Isotropy hypothesis may represent the 
small structures quite well. Moreover, looking at the kinetic energy 
spectrum we My assume that three Zones are defined, the production zone 
(low wave numbers) , the inertial zone, and the dissipative zone (high wave 
numbers) . By this way, a time scale based upon the time that the gross 
structures transfer their energy to the dissipative structures may be defined. 
This time scale Is of the order of k/e. So, the processes are simplified, 
and the e equation may be written In the transport form where the sources 
and dissipation texms are relatod to those existing In the k equation, by 
the time scale. Then following Lumley and Khajeh-Nourl (1974), the equation 19 
reduces to 


(pe)+(peu^),j + (pe^),^ - ^ 



The "diffusion term" may be approximated by a gradient formulation (Jones- 
Launder) . Then with the formulation for the eddy viscosity, the model is 
complete; the constants <7|^, a^, Ce 2 appearing in equation 20 are those 

found by Jones and Launder (1972) 




+ P - pe 


(pk)+(pku^) - g^f 

(pe)+(peOJ), - g^^ & e, ] + c f P - c 4 

J Pe ^1 ^2 k 


( 20 ) 


The formulation for the eddy viscosity proposed by Jones-Leunder is given by, 

^ z 
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This focBula assuaet s local squillbrlum Ibetwesn the production and the 
dissipation of turbulence. But In the case of a rotation, the production 
and the dissipation evolve differently becauae of the differences in the 
range of wave nusber the two processes are occurring. As the production is 
principally located in the low wave number zone, the effect of the rotation 
should be more important on the production than on the dissipation which 
can be seen as nearly laotropic. Thersfors, the relation giving the eddy 
viscosity should be corrected in the csss of a rotation. This is being 
presently carried out. 

We may remark at this point that no direct effect of rotation appears 
in the dissipation equation derived for high Reynolds number flows. That 
can be related to the weak anisotropy of the dissipation in that case. The 
analysis given above has been carried out earlier, except for its representa* 
tlon in a generalized tensor form. 

k-e Model for Low Reynolds Wuiaber Flows 

First of all, we may want to define what we call a "low Reynolds number 
flow". If we are Interested on the statistical properties of the turbulence, 
it is useful to define local parameters, such as local Reynolds number based 
on the turbulent quantities, Richardson numbers, etc. The Reynolds number 
which is interesting to characterize the turbulence behavior is based on 
the kinetic energy and its dissipation rate and is written as follows: 

^2 

R^---Re 

The fact that a flow will be of the high or low Reynolds number type will 
depend on the value of R^. 
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For low Rtynoldi fiu«b«r (lowo th« vlooouo dlffuoion My not bo nogloctodi 
Mtaovor vary naar a wall tha dlailpatton 1» not lantroplo. Many of tha 
tania wa wara naglictini nay ba Inportant vary naar a wallt Tharafora wa 
MUBt go back to aquatloni 18 and 19 and analyaa tha tarma which wara 
naglactad baCora. 

* 

Equation for tha Turbulant Klnatic Enarav k 

Tha production and turbulant dlffualon tarma ara modalad alnllar to tha 
high Raynolda nuabar and tha affact of low Raynolda numbar will appear 
through tha addy vlacoalty> Tha only tara which naada to ba aodalad la tha 
vlaeoua tarm (up^) i j . Thii tara la aqual to + U(u’^up 


Canarally tha value of tha dlaalpatlon at tha wall la diffarant froa 
saro. Following tha arguNint of Jonoa and laundar (1972) , It la uaaful to 
aaka tha dlaalpatlon e aqual to aaro at tha wall for Iniprovamant In calcula» 
tlon« It la nacaaaary, than, to know tha value of tha dlaalpatlon at tha 
wall* Wa may write tha kinetic anargy aquation at tha wall, than aquation 18 
bacomaa 


(pe) 


(p) 





J (p) 


And introducing tha following hypothaala: 

1 , ^ *** 

2. If y la the normal direction to tha wall, wa may write 

2 X 

u' w a(t)y; v’ “ b(t)y ; w* ■ c(t)y, Tharafore vary near the wall (y < 8) 

“T — ^ .2 

tha kinetic anargy la k « (a + c*"#! *y j then the flrat darlvatlva of k 
is aaro at tha wall, while the aacond darlvatlva la a conatant . 

u2 3/2 

3 , • G F p “ with e •V k ' /i and ft la a mixing length then 

THUS 

2 

y and the flrat darlvatlva of ahould be aaro, while tha aacond 
derivative ahould be a conatant. Then (pG)^pj * a W “^'^T ij^i^(p) ' 


Mor«ov«T| If ua« th« corraction function Fy for tha aublayar In tha 
ralatlon for it appaara that va laay naglact tha aacond tarm. Tharafora 
(pe) (p) “ I U '»hich ia a eonatant. 

thia can ha writtani following Jonaa-Laundar (1972), aa 


(oe) - M u .1/2 


Which can ba aaan aa a ganaraliaation to tha fomulation propoaad by Jonaa 
and Laundar for a two diawnaional boundary layar> 

1 / 2^2 


(p) 


Tha kinetic energy aquation reduces than to 

(pE)+(pRu^) + p - pe - 

It J 


10 ,« „ tl/2 cl/2 


( 21 ) 


Equation for Diaaipation Rate 

The effect of rotation may be Important in low Reynolds number flows 
(near a wall) unlike the high Reynolds number case (particularly in the case 
of turbomachines). This effect must be modeled in the equation. But this 
term is not the only one to be retained in the equation, some other terms 
which are related to the anisotropy of the dissipation at low Reynolds 
numbars are of the same Importance. 

Two kinds of terms may be discerned, which are important at low Reynolds 

numbers (see Table 1 and equation 19). 

- 1/2 

1. Terms of Order R^ 

These terms come from the interaction of the dissipation with the 
mean gradients and with the rotation velocity. They also arise from the 
terms representing the production of velocity gradients by stretching by 
fluctuating strain rate and representing the destruction of these gradients 
by viscosity on the other hand. 
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Th« tams falatad to Man quantitlan ara 

% 

Than following Liwlay (1970) thla tars can ba raducad to: 

- f u /f [» + s(l + 2 

i.J 

In tha caaa of a turboMchinary with H ■ (0, 0, n.), tha coafflclant 

p ■ 

r ^4 4« 

|l -¥ 2 -lEi— 




nay ba In fact thla tavm appeara to ba laportant 

only vary near a wall. Thacafora* in thla ragion tha boundary layar approx!'^ 
natlona ara valid and tha coafflclant nay ba approxlnatad ualng tha only 
valoclty gradlanta which ara laportant. Wa My^ than daflna a gradlant 
Rlchardaon nuabar aa followa: 

*ic - 


Hanca, 


e. .dF 

I + 2 -iEJ— - 1 - R 


U 


ic 


I.J 


If we aaaume that a relation of the gradient type exists between the 
Reynolds stress tensor and the strain tensor and that the energy spectrun 
Is not different fron an aqulllbrlun one. then we aay assume that this term 
is proportional to the production, then the terms related to mean quantities 
can be proportional to 

3 k ^ ^1 ^®T’ \c^ 
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Th« t«c»s related to fluctuating strain rate are 

- «uvg"J s;« s;, „ 

Following Lunley (1970) we approxlaate these two terns as follows: 
8*^^ jUjJ a (B + a^ v/e) 

- 1/2 

Then the part which Is of order R,j^ may be approxlnated by: 


4 f P Pj <V 


- 1/2 


So th« total contribution of tarns of order may ba raducad toi 

f ■■ "<*T- *1C> 


where the function F must be established. 


2. Terms of Order 


«-l 


There are two terns of this kind, the first one represents ''exactly" 
the diffusion processes by molecular viscosity and does not need to be 


lie 

modeled; the second term is -4u • We may note that the correlation 

between and u'^ should be weak because each term belongs to a different 
range of wave number. A first approximation then, is to consider the term 
®ik'**^ as a flux of by u*^ and to Introduce the gross hypothesis that 
this flux is related to the gradient of the mean strain and then: 


-«M sfj = 4 mVj 

However this term should be of little Importance in the major part of the 


boundary layer. 
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It Is rsBsrksbls that this tsm can ba consldarad as a ganarallsatlon 
of a tara proposad by Jonas and Laundar (1972) for tha calculation of a two 
dlsMinslonal boundary layar. Tharafora, tha dissipation aquation raducas to; 


(pe)+(pe8^>. 





R )— P 
“ic'k ^ 



“ '’l 


( 22 ) 


Tha constants Cg2 **** ** Raynolds 

nuabar casa. Tha constant reduces to tha value 2, In tha case of a two 
dlaanslonal boundary layer (sea Jonas-Launder (1972)], Than* tha daaplng 
function Fg aust be established. According to the analysis given above 
should ba a linear function of the Richardson nuiibar and an exponential 
function of the Reynolds nuaber. Than the fora of F^ should ba: 


»ic> • 1 - f<V(i - “ *ic> 


Eddy Viscosity Law 

We need a foraulatlon to relate the Reynolds stress tensor to the 
strain tensor, through an eddy viscosity coefficient. This hypothesis means 
that [Tennekes-Lumley (1972) ] the gross structures get the tendency to be 
oriented with the principal axis of the strain tensor to extract more energy 
to the mean field, and that the turbulent structures are approximately 
convacted within the mean velocity direction. Therefore, the main hypothesis 
to derive the viscosity law assumes that the turbulent viscosity is Isotropic. 
But for three dimensional flows and particularly for boundary layer flows, 
the velocity vector U and ?U are not aligned in general. Therefore, the 
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Isotropic «d4y visdosity It not adapted to pradlct tha bahavlor of tho Raynoldt 
atraaaaa. The beat way to avoid thla problem would be to solve the complete 
set of Reynolds stressas equations. However, as first step attempts will 
be made to derive an expression for eddy viscosity. In fact. In Its most 
general form the eddy viscosity It a fourth order tensor, and for the 
particular case where only two directions In the flows are of equal Importance 
(that Is the case for blade boundary layers) we may introduce two different 
eddy viscosity coefficients which can take account for the non-lsotropy of 
the flow. This has been done for a mixing length hypothesis by Kootllln and 
Lockwood (1974). Using this concept, the anisotropy Introduced by the rota- 
tion could be included in the definition of the viscosity lav. 

An attempt is presently under development to take account for the effect 
of the rotation in the eddy viscosity law. 

Conclusion 

V The analysis of the kinetic energy and the dissipation equation has 

been performed, and some remarks can be made at this time. 

Firstly, an order of magnitude analysis of equations 18 and 19 showed 
that the rotation does not affect the dissipation rate explicitly in high 
Reynolds number flows. In fact, in such flows the dissipation is nearly 
isotropic and it is logical that the rotation does not affect e because of 
the analogy of each direction. At low Reynolds number, the dissipation 
becomes non-isotropic, both the effects of Reynolds number and Richardson 
number may be important in this case. 
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Secondly, It appoers that the af facta of rotation are aora laportant 
on tha production tarma than on tha othara. In fact, if wa axanlna aqua- 
tion 17 it ia aaay to ahow that the Raynolda atraaaaa nay be greatly 
affected by tha rotation while the kinetic energy and Ita diaaipation rata 
are not ao nuch affected. Then it ia evident that the noat Important relation 
which controla the calculation la the eddy vlacoaity coefficient. 

Thirdly, we may make another remark which la alao related to the 
preceedlng covaent. The analyala of the diaaipation equation ahowa that 
the effect of rotation ahould affect the ’’production" term inatead of the 
’’diaaipation" tern in the equation for e. Moat of the modela at the preeent 
tine account for the rotation through the "diaaipation" term [aee Launder, 
Prlddln, Sharma (1977); Howard, Patankar, Bordynuik (1980)]. Nevertheleaa, 
our analyala aeema to be in accord with the remark which waa made by 
Launder et al. (1977) that the correctlona might have been better made on 
the "production" term of the e equation inatead of on the decay part. 

A computer code haa been written to check the model preaented here, in 
almple caaea auch aa flat plate boundary laycra, boundary layera on rotating 
cylinder. The program la baaed on the Patankar-Spaldlng procedure (1970) 
and la uaed to aolve the parabolic tranaport equatlona for the velocltlea, 
the kinetic energy and ita diaaipation. It la operational for two dimenalonal 
boundary layera on flat platea, and is being modified to calculate boundary 
layera on rotating cylinder. Theae calculationa are performed in order to 
check the aaaumptiona made for low Reynolda number modeling. Then the model 
could be included in the three-dlmenalonal computer codea developed at Penn 


State. 
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3. MEASURIMEMTS OF THE THREE DIMEHSIOWAL FLOW 
INSiro AN AXIAi FLOW COMPRESSOR ROTOR PASSACE 

ExptI— nf I Protr— 

Thtt MasurtaMiits rcportad In this report war* parfomad using tha Axial 
^ Flow Coaprassor facility In tha Dapartaant of Aarospac* Englnaarlng. 

Tha flow. Is survayad across tha antlra passage at six axial locations 
(ona upstraaa, four Inside tha rotor passage, and one downstraaa) and at 
five radial locations. Maasuraasnts ware taken at five radial locations 
(R « 0.58, 0.67, 0.75, 0.83, 0.918) at aach of the following axial stations 
Inside the passage: Z ■ -0.5 (upstraan) , 0.26, 0.3, 0.73, 0.97, and 1.06. 

Tha data Inside and downs trean of the rotor ware acquired with a flva- 
hola probe, rotating with the rotor. The data upstram vara acquired with 
a stationary five-hole probe. All the measurements war* taken at the design 
flow coefficient (p ■ 0.56. 

The results of the measurements are compared with tha predictions from 
the Katsanls and McNally (1977) computer program. 

Typical Results 

A small sample of the experimental data and some comparisons with the 
predictions from the Katsanls and McNally (1977) program are presented here. 

The measured and predicted blade to blade distributions of the axial 
(WZj^ tangential (NT) and radial (WR) relative velocity at R ■ 0.918 (near 
the tip) and at Z > 0.26, 0.5, 0.73, 0,97 are shown In Figures 1 through 4. 
All the components of velocity are normalized with reapect to the tip speed 
velocity. The predictions for the axial velocity are quite good at the 
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latfc two axial atationa (l.a., Z ■ 0.73, 0.93). Ac Cha first two axial 
stationa, tha agraaaant is only raasonabla and tha slopa of tha pradictad 
and naaaurad profilas ara of opposite sign. 

Tha pradictionS for tha tangantial valocity araqulta good axcapt at tha 
last axial acation Z ■ 0.97. Nayarthalaas , at this location tha shapa 
of tha pradictad and aaaaurad profilas ara strikingly idantlcal. The diffar- 
anca in nagnituda coms fron tha fact that cha coda cannot predict tha flow 
naar tha trailing adga accurataly. This is bacausa it cannot incorporate 
tha KuCta condition correctly, since it is based on an invlscid analysis. 

Tha program does not have the ability to predict blade to blade distribu- 
tions of radial valocity. The predicted radial valocity profiles that are 
shown in Figures 1 through 4 ara conputad by assuming that tha meridional 
angle calculatad in tha $2 surface solution is constant along the entire 
blade passage. 

The measured blade to blade distributions of Che relative total pressure 

and the static pressure are shown in Figures 5 through 8. Both the total 

2 

and the static pressure are normalized by 1/2 p U^. 

The static pressure profiles are fairly linear with the loading 
decreasing as we go from Z ■ 0.26 to Z - 0.97. 

The relative total pressure profiles are almost flat . The mean value 
of relative total pressure is almost constant. The difference of the mean 
relative total pressure at Z - 0.26 and Z - 0.97 is less than 1 percent. 

This indicates that the losses in the Invlscid part of the flow through the 
blade passage are very small. 

The measured profiles of the axial, tangential, and radial relative 
velocity downstream of the rotor are shown in Figures 9 through 13. The 
axial location is not the same for different radial locations due to the 
geometrical constraint on the probe. 
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Th« velocity defect le Mxiatia et the locetion Z ■ 1.049, R * 0.587 
end ■Inlaua et the locetion Z » 1.085, R ■ 0.918. 

The redlel velocity Is less then 10 percent of the tip speed velocity 
et ell locetlons end noetly outwetds. This egress veil vlth Mesureaents 
teken by Drlng et el. (1981) end Devlno end l4ikshalnereyene (1981) . Ssuill 
Inverd redlel velocity et the pressure side region esn be seen et the 
redlel locetlons neerest to the hub end ceslng while In the ald-spen region 
the redlel velocity Is everywhere outwerde. 

Inwerd redlel velocity et the pressure side region Is Induced by the 
shed vortlcity. Shed vortlclty results froa the gradient in blede loading 
across the span. Since the loedlng in the aidaspen is relatively unifora, 
there la no strong shed vortlclty in thet region so inward redlel velocity 
cennot exist. This Justifies the present experimental results. 

The aeasured profiles of the relative total pressure end the static 
pressure downstreaa of the rotor ere shown in Figures 14 through 18* As it 
was expected the relative total pressure profiles follow the behavior of 
the axial and tangential velocity profiles. 

The static pressure profiles show an Increase of the static pressure in 
the wake region. This comes in contrast to the classical assumption of 
constant static pressure across the wake. Slallar observation has been 
reported by Lakshalnsrayana and Uavlno (1980). 
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Mo— nttt« Equation In 9 Dlfction 
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APPEHDIX 2 


Llf ratuf Survty on th« Ef f«eti of Curvtuf 
and notation on Turbulanca Structura 

In thia aaction a litaratura aurvay on axparinantal and analytical 
Invaatlgatlon of tha curvatura and rotation affacta on turbulanca la carrlad 
out. 

Exparlnantal raaulta provlda Inforoatlon on the phyalcal bahavlor of 
tha turbulent flowa to enable nodalllng the effecta of the different atralna 
auch aa mean atralna, curvature effecta, rotation effecta, etc. The triple 
velocity correlatlona, preaaure-veloclty correlatlona, dlaalpatlon which 
appear In tha Raynolds atraaa aquation, are difficult to meaaure and, 
hence, the modeling la mora or leaa empirical and la baaed generally on vary 
reatrlctlva aaaumptlona. For axampla, tha turbulanca modala era of tan baaad 
on the homogeneoua fluid propertlea, local laotropy of dlaalpatlve acalea, 
high Reynolda numbera which lead to relatively almple modela. On the other 
hand, experimental work provide reaulta that can be uaed to check different 
turbulence modela. Kevertheleaa, aome of the experlmenta do not provide 
aufflclent information to check all the modela In detail. In the paat few 
yeara. It aeema that a great effort haa been made to provide experimental data 
including all the Reynolda atreaa tenaor terma. iloat of the data available 
la for aimple two dlMnalonal ahear flowa (free ahear flows, boundary layers, 
Jeta, wakea, ducta) . For the complex flows (following the definition of 
Bradshaw (1977) for complex flows) which are the most coimon In engineering 
practice, very little data and few models are available. 

In thla section, more eaphaals la given to the effect of rotation which 
la the principal effect we are interested in. 
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Curvatttf Effaete 

Ex Tparl— ntal Data 

Caatro and Bradahaw (1976) hava carrlad ouC exCanaiva ona>point naanujca- 
■anta In a highly curved nixing layar to dataxialne the effecta of aCreanllne 
cucvaeura on a ahaar layer > the principal affect la to dlalnlah the Raynolda 
atraaa tanaor components whan the angular monantum of the naan flow Inoreases 
with radius of curvature and to Increiia thaif quantities in th e opposite 
situation. The most spactsculat feature of the measurements is that the 
components of the Reynolds stress tensor after decreasing in the highly 
stabilizing curvature region, rise rapidly and overshoot the plane>layer 
value farther downstream before finally decreasing. This indicates the 
Inadequacy In current modeling for shear layers such as the use of the 
shear layer thickness to provide a length scale, the rotational Invariance 
of turbulence models based on second order transport equations, and the 
gradient diffusion hypothesis for turbulent transport. 

Changes In turbulence properties can occur even with small curvatures 
as can be seen In the experiments of Hoffmann and Bradshaw (1978) on a 
turbulent boundary layer with a mild longitudinal curvature, and of Hunt 
and Joubert (1979) in duct flow. Data are also reported by So and Mellor 
(1975) who show that Reynolds stress Increases and three-dimensional 
vortices exist over a concave surface. The experiments of Margolls (1963) 
deal with the unstable effects In a curved mixing layer. The latest paper 
by Nakano et al. (1981) covers the effects of stable and unstable freestream 
on a turbulent flow over a concave surface, where different shear flows are 
provided at the Inlet to the curved section. Three-dimensional longitudinal 
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vortiets an found. Xn an unafcabla fraa atraaa, turbulant; Intanaltlaa in 
tha boundary layer as well aa In the free atrean are Increaaed. Alcarar 
(1977) atudled the wall Jet developing on a conatant-radtua convex curved 
aurface. The flow waa nearly two-dliaenalonal, but the curvature waa large 
enough to produce Maaurabla effecta on turbulence. 

Analyala 

The effect of stabilizing or deBtabilizing forcea on the turbulent 
motion of homogeneoua flulda In flows along curved aurfacea was first dla- 
cussed by Prandtl (1930) . Early experiments gave evidence that the Intensity 
of turbulence Increases on concave surfaces and decreases on convex surfaces. 
Similar effecta occur, when tha gravity field acts on a flow of variable 
density. Gdrtler (1959) has pointed out the analogy between the two effects. 
Bradshaw (1968) proposed a formal algebraic analogy between meteorological 
parameters, such as Richardson number and the parameters describing the 
effect of curvature or rotation on turbulent flows. Seml-quantltatlve use 
of the analogy shows that the effects of curvature are appreciable if the 
shear layer thickness exceeds roughly 1/300 of the radius of curvature. The 
main result of the paper Is the Introduction of the Richardson numbers for 
streamline curvature and rotation which represent the ratio of the "buoyant" 
production to inertia production. Bradshaw also proposed a form of the Monin- 
Oboukhov formula for the change of mixing length with Richardson number for 
curved flows. The development of multi-equation models for curved flow began 
nearly 1975 with Mellor who In fact used the same stress closure approximations 
for buoyant flow and curved flow to produce a modifying function for the 
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•ddy vitcoilty. iKwln and Smith (1975) •Impllflad th« atraas-clotura of 
Laundat at al, (1975) to calciilata tha davalopmant of boundary layara and 
wall JatB on curvad ■urftcei. The noet Important riiult waa that tha 
obaat’vad curvatpra af facta could be accounted for by .the relatively muiII 
production tarma appearing In the Individual Reynolda atraaa aquatlona. 

Launder et al. (1977) proposed a two-equatlon model (k-S)> The energy 
production rate due to curvature appeara Ih exact form in the einergy equation » 
but the curvature effect la modeled empirically In the dlaalpation equation. 
Glbaon (1978) d<.veloped an algebraic Reynolds stress models following the 
idea that effects of curvature on heat transfer could only be accounted for 
by modelling the Reynolds-stress and heat->flux equations. The model Is 
derived from that developed for the buoyancy affected turbulence. The 
Influence of the wall is Introduced in the modeling of the fluctuating 
pressure. It Is shown that the effects of streamline curvature on heat 
transfer are probably significantly less than on the shear stress. That Is 
an Important result which suggests that the use of a constant turbulent 
Frandtl number in prediction methods may provide misleading estimates of 
the heat transfer from curved surfaces. The algebraic Reynolds stress model 
may be coupled with a one equation scheme (k) or two-equation (k-e) models 
to provide a length scale. Recently Townsend (1980) Introduced the rapld- 

— ."5 

dlstorsion approximation to predict the streamwlse variation of uv/q , 
considering that the complex distortion involves taking account for the 
history of the distortion. Two remarks are to be made. Firstly^ this 
approach may not be suitable for mild curvature surfaces, where the time 
scales of turbulence motion and distortion can be of the same order. 

Secondly, this method predicts only stress ratios, and If It were to form 

part of a calculation scheme, other equations would be needed to determine 

% 

the Intensities. More recently, Gibson and Rodl (1981) have proposed a full 
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Rt)molds atrsds nodal for two-dlnonalonal curvad flowa and for high Raynolda 
nunbara. Tha nodalllhg is pfincipally baaad on Laundar at al. (1975) nodal, 
ualng the alnpla gradient diffuaion hypothaaia for tha triple velocity 
correlation due to Daly and Harlow (1970) . In adapting the preaaure-atraln 
correlation to curved flow, the mean ahear production aa wall aa the extra 
atrain due to curvature are Included, though the effecta of curvature are 
Inpllcltly Introduced in the modeling. The modeled e equation. In the form 
originally propoaed by Hanjallc and Launder (1972) appeals to be quite 
adequate, only a logical change has been done In that equation. The energy 
production due to mean shear In a simple flow la replaced by the total 
(shear and curvature) production. However, recently the performance of the 
modeled equation for e has been questioned (Launder and Morse (1979), Hanjallc 
and Launder (1978)] for cooplex flows. So, one of the most updated works to 
be done should be the modeling of the exact equation. Particularly when 
curvature and rotation effects are present. We may mention the latest 
works of Gibson, Jones, Younls (1981) on curved wall boundary layer and, 

Gibson and Younls (1981) on curved wall jet more particularly on a convex 
surface [Alcaraz (1977)]. The only difference in modeling with the previous 
work of Gibson and Rodl (1981) Is In the way of how the pressure strain 
correlation Is modeled. The proposed modification for near-wall effect by 
Launder et al« (1975) Is used. All these calculations are part of a program 
motivated by the requirement for a general prediction procedure for complex 
shear flows with density stratification, rotation and streamline curvature. 


Rotation Effaces 


Expsrlnsntsl Data 

It was pointed out by Johnston Q970) » that there ara two baalc ef facta 
of rotation. If cogqionents of the eorloHs acceleration la parallel to the 
surface on which the layers are growing, secondary flows will tend to 
develop In the mean flow field of the layers. If a conponent of the corlolls 
acceleration is perpendicular to a solid surface, soae stabilising and 
destabilizing effects are observed In the turbulence structure Itself. Both 
effects are believed to be Important In the flow fields of tsentrlfugal 
Impellers. 

The experimental data available isay be classified In two categories. 

1. Free shear flows In rotating frames (e.g., wakes), and 2. Wall shear 
flows. 

In the case of free shear flows, only the experiment of Raj and 
Lakshmlnarayana (1975) is known to us. These authors give a detailed measure- 
ment of the wake characteristics behind a rotor. The results of Ravlndranath 
and Lakshmlnarayana (1980) are also available. 

In the case of wall shear flows, one suiy class the different experiments 
available in three categories ; 

1. Rotating Cylinders’ or Pipes 

2 . Centrifugal Turbomachines 

3. Axial Turbomachines 

1 . Rotatina Cylinders In Axial Flow 

Almost all the experiments designed In this case are nearly the same. 

We may siention here the works of Blssonnette and Mellor (1970) , and 
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Lohnuinn (1973) . In each of these experiments the mean properties 
and the six components of the Reynolds tensors are measured, but the major 
effects due tc nhe rotation occur very near the wall and the turbulent 
quantities are not measured In this part of the bourtdary layer. 

Arzousumlan et al. (1981) has provided data on an axially rotated cylinder 
with rypeclal emphasis In the region very close to the moving wall, neverthe- 
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less, only the following turbulent quantities are measured: u* , w' , u’w’ . 

The turbulent stresses u'v' and w’v' are derived by Integrating the corre- 
sponding momentum equations from mean velocity measurements. Nakamura et al. 
(1981), In their study of a three-dimensional turbulent boundary layer on a 
spinning thin cylinder In an axial flow, have given some results on the 
mean properties. A universal law for velocity distribution Is also derived. 
The experiment on the three dimensional boundary layer developing In an 
axially rotating pipe of Murakami et al. (1981) show that the flow Is affected 
by two counter effects. One Is a destabilizing effect due to on Increase 
of the relative velocity at the wall caused by the pipe rotation, and the 
other Is a stabilizing effect due to the suppression of turbulence by the 
cent&'lfugal force. The dominant effect depends on the Reynolds number and 
the rotational speed. 

2 , Models of Centrifugal Turbomachines 

The effects of the corlolis forces In such machines are of great 
Importance. The experiments of I-Man Moon (1964) provide some results on 
the mean quantities as well as on the following components of the Reynolds 
stress tensor: u’ , v* , u'v’ .for a rotating speed of 165 rpm and 

0.48 < R^ < 1.92. Moore (1973) carried out similar measurements but provided 
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only Man valoclty profilaa. JohnaCon at al> (1972) In thalr axparlMnt on 
fully davalopad turbulant flow In a channali which la rotating about a 
apanwlaa axla, obaacvad thraa atablllty ralated phcnoMna explained In thia 
paper. They pointed out the Rlchardaon number la an appropriate local 
atablllty paraMter. For example, local effecta of rotational atablllzatlon, 
auch aa reduction of the turbulent atreaa In wall layera, can be related to 
the local Rlchardaon number In a simple way. In thla experiment, the authora 
alao give aome quantitative data on mean propertlea. The paper of KoyaM at el. 
(1979) on "the Turbulence Structure and Three Dlnenalonallty of a Rotating 
TWO Dimenalonal Turbulent Boundary l^yer" la one of the lateat reaulta 
available. Thla paper haa not been aurveyed at thla time. We may also 
Introduce the partial reaulta of Bertogllo et al. (1980) In a centrifugal 
teating machine which give aome Information > on the flow In actual Impeller. 

3. Model of Axial Turbomachines 

For auch cases, only few reaulta are known, in fact the only available 
data are the onea provided at The Pennaylvanla State Unlveralty—>the experl> 
ment of Lakahminarayana et al. (1972) on the turbulent boundary layer on 
a rotating helical blade. However, no turbulent tneaaurementa are Included 
In these results. The more complete experiment wasAnand and Lakahminarayana 'a 
(1975, 1978) on the four bladed rotating helical channel. Some results on 
turbulence quantities show that the radial component of turbulence intensities 
is higher than the streamwise component due to the effect of rotation. More- 
over, the flow near the annulus wall is found to be highly complex. The 
turbulent shear stress measurements show that in three dimensional rotating 
turbulent boundary layers, all three correlations are of the same order of 
magnitude inside the boundary layer. A deviation Is found between thb stress 
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tensor and the strain tensor, that is one of the aost iaportsnt features in 
three-disMnsional flotrs and represents the anisotropy of the turbulence in 
such cases. 

Analysis 

Bradshaw (1969) carried out a simple analysis for both curvature effects 
and rotation effects and pointed out that the same arguments could be applied 
to shear layers In rotating fluids and curved flows. The author Introduced 
the gradient Richardson number, the flux and stress Richardson number by 
analogy with buoyancy analysis. A more detailed qualitative analysis has 
been carried out by Johnston et al. (1972) for the case of a rotating boundary 
layer. In their study, they focused attention on the production terns which 
appear in the Reynolds stress equations and which are due to the Interaction 
by the mean flow and due to corlolls effects, explicitly . In summary, this 
examination of the production terms lead to conclude that In wall layers the 
sign and the magnitude of rotation effects might be controlled by a local 
dimensionless parameter. This parameter may be related to the gradient 
Richardson number proposed by Bradshaw and Is very useful to characterize 
the stability or Instability of the flow submitted to a rotation. It appears 
that the first "theoretical" investigations and prediction procedures 
accounting for all the effects of the flow situation, namely three- 
dlmenslonallty and turbulence, are quite new. In fact the first attempts to 
calculatci chree dimensional boundary layer are based on integral methods 
such as Moore (1973), Lakshminarayana et al. (1981). The first differential 
calculations seem to be those of Majumdar et al. (1977) and Sharma (1977). 
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In tht foraar cnsc» the author* ar* aolving tha thra*«diMnalonal boundary 
layar with th* nathod of Patankar-Spaldlng (1972) , which ia a marching 
integration tachniqua for thra*odia*naional boundary layar problana. Tha 
turbulanc* modal which ia uaad in thia papar ia th* wall known k-e modal 
without any modification* for tha rotation affacta. In th* latter caaa* 
a.g. Sharma'a paper* th* calculation over an axially aplnnlng cylinder la 
performad The problem ia two«dim*naional and aolvad with the two- 
dlmanalonal marching procedure of Patankar«Spalding (1970) , Tha k>^ modal 
ia uaad to calculate tha boundary layer up to th* wall. An additional term 
la Included in the tranaport aquation for e to account for the curvature 
effect on the dlaalpatlon rate. Thia term is aeen aa a correction to the 
"dissipation" of dissipation rate e by the Richardson number based on the 
turbulent time scale k/e. Nevertheless, this kind of model cannot account 

for the non-isotropy of the stress tensor. Moreover, no correction for the 

♦ 

rotation effects are Introduced. In fact, these effects are not very 
Important for the calculation of flow on spinning cylinder. On the problem 
of the rotating cylinder, an approach similar to Sharma’s Is due to 
Spitz (1980) who solved the same momentum equations as Sharma with the same 
numerical procedure. In this work an attempt is made to Include the non- 
Isotropy of the Reynolds stress tensor by taking account for two eddy 
viscosity factors which represent the difference between the two principal 
directions. The turbulence model Is based on the mixing length hypothesis. 
We may also mention the work of Coustelx and Aupol (1979) who mad* the 
same as3unq>tlons but used a k-e model. But In each of these works no 
rotation effects were Included. 
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On« of tho firat analysla of tha rotaClon affaeta on tha turbulanca 
quantif.laa haa baan carrlad out by Raj and LakahBlnarayana (1975) and 
Raj (1975), who gava tha axact aquation In ganarallzed tanaor notatlona and 
analyzad tha turbulant procaaaaa in tha waka of tha rotor blada, but no 
ealculatlona warn carried out. Lakahalnarayana and Raynolda (1979) carrlad 
out r. qualitative analyala of tho affaeta of rotation on turbulanca In the 
naar waka of a rotor. Thla analyala Indlcataa that tha rotation haa aub* 
atantlal affaeta on tha atructura of turbulanca. Such aa radial conponant 
of Intanaltlaa la higher than tha axial and tangential In the near wake and 
decay more rapidly than the othera. The radial conponents of tha atraaaaa 
arc generally higher than thoae of a corraapondlng non-rotating caaa. 

Another analyala baaed on apactral calculua has been carrlad out by 
Bartogllo at al. and Bartogllo (1980) to atudy the effects of the rotation 
on an honogeneoua turbulent field. The effects of stabilizing and destabilizing 
due to Coriolis forces are observed, but this kind of approach does not 
account for the non-linear and Inhooogeneous terms In the cooqiutatlon. 
Nevertheless, some Important results show that the pressure-strain correla- 
tion In rotating frames may have to be modeled carefully. Another fact Is 
that the use of Isotropic functions when modeling may be Inadequate, and 
some parameters, like the direction of the force, have to be taken into 
account. 

The only conplcite calculation available, which takes Into account the 
rotation effects. Is the one performed by Howard at al. (1980). These 
authors, following the work of Majumdar et al. (1977), used a modified 
procedure based on the partlally-parabollc method of Majumdar and Spalding 
(1977) . They solve the three-dimensional boundary layer equations within 
a rotating frame. The turbulence model is based on the k-e model. Three 
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•odlflcationa to the basic k-e modal ara Caatad which includa affacts of 
Coriolia forca on tha turbulant anargy and dissipation rata. Tha first modal 
is basad on the work of Wilcox and Chaad>crs (1977) and taka account of the 
Coriolis affacts in both tha k and e aquatiohs assuming that the kinetic 
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energy is proportional to u* . Tha second and third models ara basad on 
a study of curved boundary layara by Launder at al. (1977) » where the k 
aquation is maintained in its usual form and where the Coriolis effect is 
Introduced through the e equation. Two forms of the Richardson nuad>er ara 
considered. Their conclusion is that the Wllcox-Chambers model give the 
most satisfactory prediction, and they point out the need for inclusion of 
a Coriolis modal for turbulence modification. Finally, the approach of 
Hah and Lakshmlnarayana (1980) and Hah (1981) who introduce an algebraic 
Reynolds stress model coupled with a k-e model which, though the convection 
and diffusion are not Included, take account forrotatlon and curvature 
effects . 

Comments 

The review of the different paper point out three main features: 

1. Only feW' calculations are available for the calculation of the 
three-dimensional boundary layer in rotating frames. 

2. No complete Reynolds stress model is available for rotating turbulent 
flows. Very few attempts have been made to account for the rotation effects 
in the k-e model, but none are really based on a logical analysis. In fact, 
the most up to date results in modeling are those of Raj (1975) and Hah and 
Lakshmlnarayana (1980) . The major effort, then, should be given to the 
analysis of the dissipation rate equation and to the Reynolds stress equations. 
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3, D«tall«d MaaurcMntB providing inforaation on tha af facta of 
both tha Roaaby and tha Richardaon maibar on turbulanca ara not vary 
nunaroua. So, avary naw raault would ba of graat intavaati particularly 
ii tha rotation affact can ba laolatad fron tha othar af facta. 
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Figure 1. Blade to blade distribution of axial (WZ), tangential 
(WT), and radial (WR) velocities at Z » 0.26, R - 0.9 
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Figure 2. Blade to blade distribution of axial (WZ) , tangential 
(WT), and radial (WR) velocities at Z ■ 0.5, R « 0.918 
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Figure 12. Axial (WZ), tangential (WT), and radial (WR) velocity 
profiles at Z • 1.077, R ■ 0.832 
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Figure 13. Axial (WZ) , tangential (VJT) , and radial (WR) velocity 
profiles at Z * 1.085, R ■ 0.918 
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